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We present a detailed discussion of some features of quantum mechanical metastability. We 
analyze the nature of decaying (quasistationary) states and the regime of validity of the exponencial 
law, as well as decays at finite temperature. We resort to very simple systems and elementary 
techniques to emphasize subtle aspects of the problem. 



I. INTRODUCTION 

Complex-energy eigenf unctions made their debut in 
Quantum Mechanics through the hands of Gamow, in 
the theory of alpha-decay.^ Gamow imposed an "outgo- 
ing wave boundary condition" on the solutions of the 
Schrodinger equation for an alpha-particle trapped in the 
nucleus. Since there is only an outgoing flux of alpha- 
particles, the wave function 'ijj{r, t) must behave far from 
the nucleus as (for simplicity, we consider an s-wave)^ 



-iEt-i-ikr 



(r oo). 



(1) 



This boundary condition, together with the requirement 
of finiteness of the wave function at the origin, gives rise 
to a quantization condition on the values of k (and, there- 
fore, on the values of i? = It turns out that such 
values are complex: 



En 



It follows that 



-r„t+if„r 



jr„/2. 



(r — > oo). 



(2) 



(3) 



Thus, if r„ > 0, the probability of finding the alpha- 
particle in the nucleus decays exponentially in time. The 
lifetime of the nucleus would be given by t„ = 1/rn, and 
the energy of the emitted alpha-particle by e„. 

Although very natural, this interpretation suffers from 
some difficulties. How can the energy, which is an ob- 
servable quantity, be complex? In other words, how can 
the Hamiltonian, which is a Hcrmitean operator, have 
complex eigenvalues? Also, the eigenfunctions are not 
normalizable, since r„ positive implies Kn positive and, 
therefore, according to (3), i)P diverges exponen- 

tially with r. 

In spite of such problems (which, in fact, are closely 
related), it is a fact of life that alpha-decay, as well as 
other types of decay, does obey an exponential decay 
law, with a rate close to that obtained using Gamow's 



method. Why this method works is the question we try 
to answer in this paper in a very elementary way.^ Thus, 
in Section II, we show Gamow's method in action for a 
very simple potential. Some of the results obtained there 
are used in Section III, where we study the time evolu- 
tion of a wave packet initially confined in the potential 
well defined in Section II. This is done with the help of 
the propagator built with normalizable* eigenfunctions, 
associated to real eigenenergies. As a bonus, we show 
that the exponential decay law is not valid for very small 
times or for very large times. This is the content of Sec- 
tion IV, where the region of validity of the exponential 
decay law is roughly delimited. 

Another topic we address in this paper is decay at fi- 
nite temperature. This is done in Section V, where we 
study, with the help of an exactly solvable toy model, 
how the decay of a metastable system is modified when 
it is coupled to a heat bath. It is shown that, under suit- 
able conditions, the decay is exponential, with a decay 
rate T given by the thermal average of the r„'s. 



r 



-,3e„ 



(4) 



Although this result appears to be rather obvious, in fact 
it is not: the decay of a metastable system is an intrinsi- 
cally non-equilibrium process and, so, there is no a priori 
reason for the decay rate to be given by (4). Finally, in 
Section VI, we discuss the concept of "free energy of a 
metastable phase," a point where we think there is some 
confusion in the literature. 

The results and ideas presented here are not really 
new, but discussions on these matters usually involve 
the use of sophisticated mathematical techniques, such 
as functional^'^ or complex analysis. ''""'^^ For this reason, 
we have tried to make the presentation as clear as possi- 
ble by resorting to very simple systems and elementary 
techniques — in fact, techniques that can be found in 
any standard Quantum Mechanics textbook. 
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II. DECAYING STATES 

In order to exhibit Gamow's method in action, we shall 
study the escape of a particle from the potential well 
given by: 



for a; > 0, 
for a; < 0. 



(5) 



Motion in the region x < is forbidden because of the 
infinite wall at the origin. The positive dimensionless 
constant A is a measure of the "opacity" of the barrier at 
X = a; in the limit A — > oo, the barrier becomes impene- 
trable, and the energy levels inside the well are quantized. 
If A is finite, but very large, a particle is no more con- 
fined to the well, but it usually stays there for a long 
time before it escapes. If A is not so large, the particle 
can easily tunnel through the barrier, and quickly escape 
from the potential well. Metastability, therefore, can only 
be achieved if the barrier is very opaque, i.e., A is very 
large. For this reason, we shall assume this to be the case 
in what follows and, whenever possible, wc shall retain 
only the first no n- trivial term in a 1/A expansion. 

To find out how fast the particle escapes from the po- 
tential well, we must solve the Schrodinger equation 

d 1 A 

i^ip{x,t) = -- il;{x,t) + - S{x - a)ip{x,t). (6) 



dt 



2 dx^ 



a 



'tp{x, t) = cxp{~iEt) ip{x) is a particular solution of this 
equation, provided (p{x) satisfies the time-independent 
Schrodinger equation 



2 dx"^ 



f{x) H — 5{x — a) ip{x) = E ip{x). (7) 



Denoting the regions < a; < a and x > a by the in- 
dices 1 and 2, respectively, the corresponding wave func- 
tions ipj{x) {j = 1, 2) satisfy the free-particle Schrodinger 
equation: 



2 ^^i(^) =•^<^i(2')• 



(8) 



Since the wall at the origin is impenetrable, </5i(0) must 
be zero; the solution of Eq. (8) which obeys this bound- 
ary condition is 



(9) 



ipi{x) = A sinkx {k = v2E). 



To determine (p2{x), we follow Gamow's reasoning^'^'^^ 
and require (p2{x) to be an outgoing wave. Therefore, we 
select, from the admissible solutions of Eq. (8), 



(10) 



The wave function must be continuous at a; = a, so that 

ipi{a) = ip2{a), or 



^ = e sinAia. 
A 



(11) 



On the other hand, the derivative of the wave function 
has a discontinuity at a; = a, which can be determined 
by integrating both sides of (7) from a — e io a + e, with 
£^0+: 



(/92(a) -vi (a) = ^<y92(a), 



(12) 



from which there follows another relation between A and 
B: 



B 
A 



-ika 



COS ka 



2\/ka 



(13) 



Combining (11) and (13), we obtain a quantization con- 
dition for k: 



cotan.a = i-^. 

ka 



(14) 



The roots of Eq. (14) are complex; when A :» 1, those 
which are closest to the origin are given by® 



kn 



2mrX 



H2a) 



.H-2A 

The corresponding eigenenergies are 



(nTT < A) 



(15) 



E^ = kl/2 « 1 ( 



= e„ 



2 \~) 
irj2. 



(2A^ 



(16) 



(Note that Kn <C k„ and r„ ^ e„; these results will be 
used later.) The imaginary part of En gives rise to an 
exponential decay of \tpnix,t)\'^ , with lifetime equal to 



Tn = l/r„ 



2(Aa)' 



(17) 



Since the corresponding value of B/A is very small (~ 
n/A), one may be tempted to say that the probability of 
finding the particle outside the well is negligible in com- 
parison with the probability of finding the particle inside 
the well. Normalizing ipn in such a way that the latter 
equals one when t = 0, the probability of finding the 
particle inside the well at time t, if it were in the n-th 
decaying state at t = 0, would be 



Pn{t) 



\ipn{x,t)\'^ dx = exp{-Tnt). 



(18) 



The trouble with this interpretation is that ImA:„ = 
—Kn/2 < 0, and so 'ipn{x,t) diverges exponentially as 
a; — > oo, since, according to (10), 

i)|^ = IBnf exp(-r„t + Knx) (19) 

outside the well. Because of this "exponential catastro- 
phe" , the decaying states are nonnormalizible and, there- 
fore, cannot be accepted as legitimate solutions of the 
Schrodinger equation (although one can find in the liter- 
ature the assertion that they are "rigorous" solutions of 
the time-dependent Schrodinger equation^^). 



2 



III. TIME EVOLUTION OF A WAVE PACKET 

We now return to Eq. (8) and write, for the solutfon 
in region 2, instead of (10), the sum of an outgoing plus 
an incoming wave: 



-ikx 



+ Be 



ikx 



Continuity of the wave function at a; = a implies 



A sin A;a = e 



-ika 



+ Be 



\ka 



(20) 



(21) 



As before, the derivative of the wave function has a dis- 
continuity at X = a, given by Eq. (12), from which it 
follows, instead of (13), 



kA cos ka ■ 



2A ., 

\-ik 

a 



-ika 



a 



Solving (21) and (22) for A and B, we find 



A 
B 



2ika 



{ka + A sin 2ka) +i2X sin ka 
{ka + X sin 2ka) — i2X sin^ ka 
{ka + X sin 2ka) + i2X sin^ ka 



(22) 

(23a) 
(23b) 



These expressions show a couple of interesting features: 

(1) \B\ = 1 for real values of k, implying a zero net 

flux of probability through x = a; therefore, milikc the 
solution found in the previous section, there is no loss or 
accumulation of probability in the well region. 

(2) 1^1 <C 1 if fcffl <C A, except if k is close to a pole of 
^(A;). in which case \A\ may become very large. 

To find the poles of A we must solve the equation 
A{k)~^ = 0, which, after some algebraic manipulations, 
reads 



cotan ka = i — 



2A 

ka 



(24) 



This is the same as Eq. (14)! Is this a coincidence? In 

fact, no. According to (23), A and B have the same poles; 
in a sufficiently small vicinity of a pole, |^| and \B\ are 
very large, and so Eqs. (21) and (22) become equivalent 
to Eqs. (11) and (13), respectively. 

Suppose that at t = the particle is known to be in 
the region x < a with probability 1; in other words, its 
wave function tp{x, 0) is zero outside the well. 



'tjj{x, 0) = for a; > a. 



(25) 



The wave function at a later time t is given by 

i;{x,t)= G{x,x';t)ilj{x',Q)dx', (26) 
Jo 

where the propagator, G{x,x';t), can be written as 



/•CO 

G{x,x';t)= / e-'''^*/'^ip*{k,x)ip{k,x')dk. (27) 
Jo 

The function (p{k, x) is the solution of Eq. (7) correspond- 
ing to the energy E = k'^/2: 



(p{k,x) 



J A{k) sin kx 



{x < a) 



(28) 



With this normalization, the (fi{k, x) satisfy the complete- 
ness relation^^ 



ip* (fc, x) ip{k, x') dk = 6{x — x'). 



(29) 



Since, by hypothesis, 'tp{x,0) = fov x > a, (26)-(28) 
give, for x < a, 



il){x,t) 



1 

2^.70 

dx't{;{x',0) sinkx'. 



For k close to a pole k 
approximated by 

A{k) 



dke-''''*/^\A{k)f sinkx 

(30) 

Kn — iKn/2, A{k) can be 

2ikr,a 



a(l + 2Ae2f=-") {k - kn) 

' {k- Kn)+ iKn/2' 



(31) 



As we have seen, |A(A:)p <C 1 if fca <C A, except at 
the resonances, where (31) may be used. On the other 
hand, if tp{x, 0) is sufficiently smooth, in the sense that 
dxip{x,Q) sinfca; — > sufficiently fast when — > oo 
(this condition will be made more precise later), then 
most of the contribution to the integral (30) comes from 
the region ka <C A. Therefore, we may approximate (30) 
by 



dke 



n 

X / dx' il){x' ,{)) 
Jo 



.ik^t/2 (/tn/A)^ sinfca: 
{k-Knr + K^/A 



sin kx'. 



(32) 



where /„ is the interval [(«„ + k„_i)/2, (k„ + Kn+i)/2] 
{n = 1, 2, . . . ; /to = 0). Because of the arguments pre- 
ceding (32), only the first few terms of the sum give a 
significant contribution to the integral. Note also that, 
since the resonance in |^(/;;)p around has a width of 
the order of Kn, and 



KnX, Knx' < Kntt w {mrf /2X'^ «; 1, 



(33) 



we can substitute k for «;„ in sin kx and sin kx' in the 
integrand of (32). On the other hand, this is not allowed 
for e"**^ since the time t is not bounded. 
Let us examine the integrals 
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{k-nnY+Kl/A- 



(34) 



If Kn <C — Kn-1, we Can safely extend the interval 
of integration to the whole real axis, and carrying out the 
integration we find 



J — ( 



1 



(fc - K„)2 + Kl/A 



1 /-OO 

V27rii J-oo J-i 



27r 



.0 
^ — c 



^^gi{72t+i(K„-iX„/2){^ 
^^giCV2t+i(K„+iif„/2)|_ 



(35a) 
(35b) 
(35c) 



Except for a region of width ^ \ft around the point 
^ = —Knt, where the phases of the exponentials arc sta- 
tionary, the oscillations of the integrands tend to can- 
cel out, giving a very small contribution to the integrals 
above. If Knt ^ \ft, such a region is well inside the 
negative real axis, therefore the second integral can bo 
neglected in comparison to the first. For the same rea- 
son, we can extend the interval of integration of the first 
integral to the whole real axis, thus obtaining 



27r 



if /2t+i{K„-iK„/2)i 



- ^ 



(36) 



(37) 



Substituting this result in Eq. (32), we find 

i^{x, t)K^Cn e"'^"* Vn{x), 
n 

where <fin{x) and c„ are defined as 



[^{Kn/Xf . [2 . rnrx 

Vnyx) ^ \ - — smK^x « W - sm , (38a) 

V I Kn yd a 



Cn= dxtp{x,0) iPn{x). 

Jo 



(38b) 



Eq. (37) is formally identical to the well known expansion 

of the wave function in energy eigenfunctions, except for 
the fact that: (1) it is an approximate result and, as such, 
subject to some restrictions, and (2) the energies En are 
complex, as a result of which the probability P{t) to find 
the particle inside the potential well at time t decreases 
in time: 



P{t) = r|^(a;,i)|2dx«^|c„pe- 



(39) 



Now wc can be more precise on the smoothness of ip{x, 0); 
roughly speaking, the smaller the value of n beyond 
which c„ = 0, the better the results above will be. 



Eq. (37) is valid only inside the potential well. To find 
the wave function outside the well, we must return to 
Eqs. (26)-(28) and make x > a: 

^P{x,t) = — / dke-''''^/^ [e''''' +B*{k)e-''''']A{k) 
27r Jo 

X / dx'il){x',0) sinkx'. (40) 

^0 



Since most of the contribution to the integral in k comes 
from the resonances, we may approximate A{k) by (31) 
and make an analogous approximation for B*{k). Thus, 
(40) becomes 



ll>{x,t) K.^dnJn{x,t); 



(41a) 



dn = -^—^ I dx' tp{x' ,0) sin Knx' , (41b) 
27rA Jo 



27rA 

Jn{x,t)= I dke-'"^'"^ 



In 



k- Kn + iKn/2 



— c.c. 



(41c) 



where /„ has the same meaning as in Eq. (32), and c.c. 
denotes complex conjugate. 

Let us concentrate our attention on the integrals 

J^n{x,t). Extending the interval of integration to the 
whole real axis, and using the same trick as in Eq. (35a), 
we find 



^n (^) ^) 



27rz 



•/ — OO 
J X 



(42a) 
(42b) 
(42c) 



As in the case of (t) , the second integral is negligible 
in comparison to the first if n^t S> Vt, or t ^ l/e„. On 
the other hand, the first integral may be approximated 

by 

/OO 
^^gi€V2t+ifen(l+x) ^ V2^e-<*/2+^'="^ (43) 
-C30 

only if Knt — X ^ \fi. 

Returning to Eq. (41), we finally obtain 



V'(a;, t) w — 27r« dn e ' 



(44) 



We see, therefore, that outside the well the wavefunction 
behaves as a superposition of outgoing waves, in the way 
postulated by Gamow. However, the exponential catas- 
trophe does not occur here, for Eq. (44) is valid only 
under the assumption that Knt — x <^ \ft. 
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IV. BREAKDOWN OF EXPONENTIAL DECAY 



Thus, Eq. (37) gives 



The evolution of the wave function requires some 

time'^^ to reach the regime of exponential decay: typi- 
cally, a time corresponding to many oscillations inside the 
potential well (i.e., t l/e„). To be more precise, even 
if this condition is satisfied, the decay is not strictly ex- 
ponential, but a sum of exponential decays, one for each 
resonance [Eq. (39)]. However, since the lifetime t„ is, in 
general, a rapidly decreasing function of n (t„ w n/n^ 
in our example), the decay becomes a pure exponential 
one after a time of the order of ti . 

On the other hand, the exponential decay does not 
last forever. After some sufficiently long time, it obeys 
a power law.'''"-'^^'-'^^ To see this, note that for t oc, 
the integral (30) is dominated by small values of k. One 
finds, then, for x < a, 



3en-r„t 



(50) 



and, as already discussed in Sec. IV, after a time of the 
order n = l/Fi, the decay would be dominated by the 
decay of the "false vacuum" — the lowest lying reso- 
nance. It follows that the decay rate is almost insensitive 
to the temperature. 

Is this conclusion correct? 

In the literature-'^^ one finds the statement that, in such 
a situation, the probability P{t) decays as 



p(t) = e-<^>*, (r) = i5^r„, 



-0e„ 



(51) 



ip{x,t) 



2w 

X 



g-ife"f/2 ^2 



A perfectly sensible question is: why does P{t) decay this 

way, and not as 



A2 



/ dx' ^l>{x',0)x' / dk 

V..>(x',0).'.° (45) Pit) = e-^^^^\ M = |E-"^-^-? (52) 

Jo " 



Therefore, the probability of finding the particle inside 
the potential well behaves asymptotically as 



Pit) 



dx'ip{x',0) x' 



2 6 



(46) 



In order to answer these questions, we shall study a toy 
model: a two-level metastable system coupled to a heat 

bath. (For instance, imagine a situation in which only the 
two lowest resonances of a potential well are excited.) 



These levels have "complex energies" Ej 



irj/2 



Comparing (46) with (39), one finds that they become 
comparable in magnitude when 



(j = 1, 2). Let ni{t) and n2{t) be the populations at time 
t of levels 1 and 2, respectively. A reasonable dynamics 
is given by the following set of equations (n = dn/dt, 
E = e2- ei): 



or, since A 1, 



t 



n 



10 lnA-h3 InlnA. 



(47) 



(48) 



ill = — Ti ni +r {n2 — e rii), 
na = -Fa na - T (ria - e-f^ m). 



(53) 



Thus, when the decay begins to obey a power law^^ 
(~ t~^), the probability that the particle is still inside 
the potential well is so small (< A~^°), that it would be 
very difficult to observe deviations from the exponential 
decay. 



V. DECAY AT FINITE TEMPERATURE 

In general, the initial state of the particle, ip{x,0), is 
not precisely known. Such a knowledge is required in or- 
der to determine the coeficients Cn in Eq. (37). Let us 
imagine, however, that the system is in contact with a 
heat bath at temperature T. Then, it is reasonable to 
assume that (/? = l/ksT, ks = Boltzmann constant) 



The first term on the r.h.s. of (53) describes the "natural" 
decay of the levels. The second term is due to the cou- 
pling to the the heat bath; it drives the system towards 
thermal equilibrium with the bath. 

Since (53) is a set of differential linear equations, solu- 
tions can be found in the form 



ni{t) = ae 



-xt 



n2(t) =be 



-At 



(54) 



-/3e„ 



-I3er, 



(49) 



The decay rate. A, must satisfy the characteristic equa- 
tion 

A^ - [Li + La + (1 + e-'5^)r] A 

+ rir2 + (ri + e-^^r2)r = o. (55) 

Although this equation can be solved exactly, the exact 
solution is not very illuminating. Instead, we shall con- 
sider two limiting situations. 

(1) r Fi, F2 (overdamping): in this case, we may 
approximate Eq. (55) by 

A^ - (1 + e-''^)FA + (Fi + e-^^ F2)F = 0, (56) 
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whose solutions are 



A±--(i + e )r|i±Wi- (i + ,-^^;)2r 



(57) 



Expanding the square root, we find (since Fi, r2 <C F) 
A+ « (1 + e-^^)T, (58a) 

(58b) 



A- 



1 + e-/'^ 
The corresponding eigenvectors are 

1 



a+ 



b- 



(59) 



If nj(0) = rijo (i = 1,2), then, at time t, we have 

'ni(t) \ e-'""nio-n2o / 1 



n2(i) 



+ 



1 + e-/5s 

«10 + "■20 



-A+t 



1 1 p--^-* 
-/3E I e 



(60) 



Therefore, even if the system is initially not in thermal 
equilibrium with the heat bath, it thermalizcs in a time 
of the order 1/A_|_ <C 1/A_. In other words, after a tran- 
sient time of the order of 1/A+, both levels decay at the 
same rate, A_, equal to the thermal average of Fi and 
F2, and their populations soon approach the Boltzmann 
distribution. 



n2{t) 
ni{t) 



-f)E 



(61) 



(2) F <C Fi, F2 (underdamping): now, Eq. (55) may 
be approximated by 



A^ - (Fi + F2)A + F1F2 = 0, 



(62) 



whose solutions are \j = Tj {j = 1,2). The correspond- 
ing eigenvectors are" 



1 




The solution of (53) is, therefore, 

n,(t) «n,oe-^^* (j = l,2). 



(63) 



(64) 



In this case, the levels decouple from each other, and 
each one of them decays with its own decay rate. The 
coupling to the heat bath is so weak that the system is 
effectively insulated. 

Now we are in position to answer the questions posed 
in this section: 

(A) If the system is ovcrdampcd (i.e., F ^ Fi, F2), 
it decays as indicated in (51). The analysis of the over- 
damped case also explains why the decay rate is given 



by (51), instead of (52). (In some sense, that is why 
one can observe a phenomenon like the Stark effect. In 
the presence of a constant external electric field, the po- 
tential which binds an electron to an atom becomes un- 
bounded from below, and so the atomic energy levels 
become metastable, with very large lifetimes if the ex- 
ternal electric field is small compared to the field of the 
nucleus. The role of the heat bath is played here by the 
quantized electromagnetic field; it is the coupling of the 
atom to it that causes the excited states of the atom, oth- 
erwise stationary, to decay to the ground state. Even for 
a highly excited atom, for which the natural lifetime is 
relatively high (some milliseconds), the situation is well 
described by saying that the atom is overdamped, since 
the ionization (decay) rate is very small if the external 
electric field is small compared to the intratomic field.) 

(B) On the other hand, if the system is underdamped 
(F <C Fi, r2), although the decay is still described by Eq. 
(37), Eq. (50) is possibly not valid. The reason is that, 
under such a condition, the system does not "know" the 
temperature of the heat bath. It would have decayed be- 
fore it could thermalize. In practice, as already argued, 
the excited states would depopulate much sooner than 
the "ground state," and so the decay rate would be very 
insensitive to the temperature of the heat bath. 



VI. THERMODYNAMICS OF METASTABLE 
SYSTEMS: A BRIEF DIGRESSION 

Finally, we would like to make a brief digression on a 

point where we think there is some confusion in the lit- 
erature. It concerns the thermodynamics of metastable 
systems. As an example of such a system, let us con- 
sider a particle interacting with the potential defined in 
Section II. Its partition function is defined as 

^ = Tre-''^ 

= / dke-l^^^'^ / dx<p*{k,x)ip{k,x), (65) 
Jo Jo 

or, alternatively, as 



/ PiE) 
Jo 



e-^^dE. 



(66) 



Clearly, the spectral density p{E) is given, in this case, 
by4 



p{E) = 5(0) 



dfc _ 5(0)_ 
dE ^ 



2E 



(67) 



We conclude, therefore, that there is no sign of metasta- 
bility in the partition function. Now, let us try to define 
a "partition function inside the well;" since this is not a 
fundamental concept, more than one definition is possi- 
ble. One such definition is inspired by (65); restricting 
the ^-integration to the interval [0, a], we have 
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Zi= dke-^^''''^ dxip*{k,x)ip{k,x). (68) 
Jo Jo 

If we make the same approximations we made in Sec- 
tion III, we can reexpress Zi as in (66), but now with a 
spectral density given by 

This kind of spectral density^^'^^ contains some dynami- 
cal information — resonant levels and decay rates. If the 
latter are small enough, the Lorentzians in (69) can be 
approximated by delta functions, and we obtain, there- 
fore, 

Z^^J2e-^'- = e-0^K (70) 

n 

If the coupling with the heat bath is strong (in the sense 
of Section V), Fi can be interpreted^^'^^ as the free en- 
ergy of the mctastablc phase. One should not cionfusc^"'^ 
this "free energy" with the true equilibrium free energy, 
jr= -(1//?) InZ. 

Another possible definition of the "partition function 
inside the well" is inspired in the complex eigenenergies 
of Gamow's method: 

n 

This definition, although almost identical to (70), pre- 
sents a new and interesting feature: the free energy F2 of 
the metastable phase is complex! Its real part is essen- 
tially equal to -Fi, and so has thcrmodynamical content, 
but its imaginary part provides dynamical information: 
if the r„'s are small, it is easy to show that 

Imi^2«-^(r), (72) 

where (F) is the thermal average of the r„/s, as defined 
in (51). That is the reason why definition (71) is so pop- 
ular in the literature, although it is no more fundamental 
than definition (68). 

VII. CONCLUSION 

In Section III we showed that decaying states, although 
plagued by the exponential catastrophe, give a fairly 
good description of the decay of a metastable state, pro- 
vided some conditions are satisfied. In fact, the main 
result of this paper is that one can compute the decay 
rate solving the time independent Schrodinger equation 
subject to the "outgoing wave boundary condition," Eq. 
(10). This is far from being a trivial result, since the 
corresponding eigenstates are nonphysical. The "effec- 
tiveness" of the decaying states in describing the decay 



may be understood by noticing^ that they are good ap- 
proximate solutions to the time-dependent Schrodinger 
equation, although nonuniform ones (i.e., they are not 
valid in the entire range of values of t and x). 

In Section V we examined another "well known" re- 
sult, that the decay rate of a metastable system in con- 
tact with a heat bath is given by Eq. (51). Although we 
have used a very simple toy model to discuss this point, 
we believe it contains the essential physics of the phe- 
nomenon, at least in the two limiting cases we studied in 
some detail. The important lesson to be learned here is 
that Eq. (51) is correct (at least in first approximation), 
provided the condition of "overdamping" is satisfied. ^^'^'^ 
At low temperatures, where only the lowest lying decay- 
ing states take part in the process, it is an easy matter to 
verify if it is so. However, as the temperature increases, 
decaying states of higher energy are excited and begin 
to play an increasingly important role in the overall de- 
cay. Since the decay rates r„ become larger with n, the 
overdamping condition eventually fails to be satisfied by 
states actively involved in the process of decay. Thus, 
one should expect deviations from the decay rate given 
by Eq. (51). Another source of deviations, not taken into 
account in our toy model, is a possible renormalization 
of the complex energies caused by the interaction of 
the system with the heat bath. This may affect Eq. (51) 
even at low temperatures, for obvious reasons. 
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